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We use Daubechies' orthonormal compact wavelets as a variational basis for the XY model in two and three di- 
mensions. Assuming that the fluctuations of the wavelet coefficients are Gaussian and uncorrelated, minimization 
of the free energy yields the fluctuation strength of wavelet coefficients at different scales, from which observables 
can be computed. This model is able to describe the low-temperature phase and makes a prediction about the 
phase transition temperature. 



1. Introduction 

The success of wavelets in analyzing complex 
signals has prompted speculation about their pos- 
sible applications to field theories and other lat- 
tice systems. We report here on the first part 
of an investigation into possible wavelet tools 
that can be applied to lattice field theories. The 
wavelet representation of a lattice field is analyzed 
by a variational principle to obtain an approxi- 
mate description of a phase transition. Wavelets 
represent nonlocal fluctuations that eventually in- 
duce the phase transition. 

1.1. The wavelet transform 

A wavelet [1,2] ip £ L2(R) is a function whose 
binary dilatations and dyadic translations gener- 
ate a Riesz basis of L^(R) so that any function 
/ G L^(R) can be expanded into a wavelet series 

oo 

/w= E E (1) 

n = — oo x'££" 

with basis functions xp'-"\x') E L^(R), 

^("\x'){x) = 2-"/'^^{2-"{x-x')) , (2) 

labeled by scale n G Z and position x E jC" on 
the n-th sublattice C" = 2"Z. 

Compact orthonormal wavelets as constructed 
by Daubechies [3] are not given analytically but 
by a powerful numerical procedure [4,5]. The al- 
gorithm acts on vectors (i.e. functions discretized 
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on a lattice), seperating the information con- 
tained in a vector into a smoothed vector on a 
coarser lattice and a complementary detail vec- 
tor. Repeated application of the algorithm strips 
the initial vector scale by scale of its detail infor- 
mation (= wavelet coefficients). Being a linear or- 
thogonal transformation, it defines an orthonor- 
mal basis consisting of wavelets (associated to the 
detail information) and scaling functions (associ- 
ated to the smoothed information). The discrete 
wavelets defined in this way approach the contin- 
uum wavelet. 

Wavelets in more than one dimension can be 
constructed by direct products. They then carry 
an additional label t = 1, . . . , 2^ — 1 that specifies 
their composition. 

2. The model 

Our model is the XY spin model defined by 
the Hamiltonian 

D 

Ti = ^ ^ {1 — cos {a{x + e^) — a{x))} (3) 

on a lattice C = in D dimensions. The field 
—IT < a[x) < TT gives the angle of the spin at 
site X, and the unit vector in direction In 
two dimension, this model undergoes a Kosterlitz- 
Thouless phase transition between a disordered 
high-temperature phase and a spin-wave domi- 
nated low-temperature phase. 
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2.1. Variational principle 

We employ a variational principle to obtain an 
approximate description of the probability distri- 
bution in the partition sum. To this end, the field 
a[x) is expanded in a wavelet basis. We choose a 
Gaussian trial probability distribution character- 
ized by the correlation function 



the internal energy 



("0 



(4) 



i.e., fiuctuations of the wavelet coefficients are as- 
sumed to be uncorrelated. The values of the vari- 
ational parameters A["\ describing the strength 
of fiuctuations at scale n, are obtained by calcu- 
lating the internal energy U (expectation value 
of 7i in this ensemble) and minimizing the free 
energy 



(5) 



Due to the Gaussian ansatz, both the entropy S 
and the internal energy U can be easily calcu- 
lated. However, the field a[x) can now assume 
any number although the Gaussian ansatz does 
not refiect the periodicity of the Hamiltonian. 
This will limit the applicability of the variational 
method to the low-temperature phase where pe- 
riodicity is not important. 

2.2. Calculation of the internal energy 

We denote the wavelet transform of the field 
a{x) by a'-"\x'), x' E C" , defined on the sub- 
lattices = (2"Z)^. To obtain the Hamilto- 
nian in terms of the wavelet-transformed field, we 
must express the difference of two adjacent spin 
angles in wavelet coefficients. This is done using 
the characteristic function of the link, 



L{x,fi){y) = 6y^^+e^ - 6y 



(6) 



and the scalar product in the space of lattice 
fields. In wavelet space, the difference is a scalar 
product of the wavelet transform l["\x, lJ.){x') of 
the link with the wavelet coefficients of the field. 
Thus the simple difference of two adjacent lattice 
sites is replaced by a more complicated, albeit 
still linear, expression, the cosine of which can be 
evaluated in a Gaussian ensemble. One finds for 



l-exp 



The quantity 



(7) 



(8) 



characterizes the power spectrum of the link 
(x, /J,), i.e., how much of the original square norm 
("power") of the link's characteristic function is 
on scale n. It can be computed numerically and 
determines how much fiuctuations at this scale 
contribute to the internal energy at this link. 

3. Results 

3.1. Free fields 

The same procedure as above can be performed 
for a Gaussian model (free field theory), corre- 
sponding to expanding the cosine and retaining 
only the first term. In this case, the minimiza- 
tion can be performed analytically yielding 



A 



(n) 



(9) 



(#„ = zf^ of sites on £"). Using the scaling prop- 
erties of wavelets, these obey a simple scaling law 
from which the correlation function C'{d) can be 
calculated: 



An 

^ 2/3 



C{d) ~ 



(10) 



exhibiting the expected scaling behavior for a 
free field theory. Note that this critical correla- 
tor results from local (noncritical) correlations in 
wavelet space. This refiects a general hope that 
wavelets may provide an approximately diagonal 
basis for some systems. 



3.2. XY model 

We have minimized numerically the free energy 
with respect to the fiuctuation strengths A^-j^\ 
At low temperatures, the model approaches the 
Gaussian model whose results serve as a start- 
ing and reference point. Spins are approximately 
aligned, and spin wave fiuctuations rarely make a 
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single spin turn by a large angle. The periodicity 
of the Hamiltonian is thus not felt by the model. 

At higher temperature, the probability distri- 
bution expands and starts to feel the periodic- 
ity. At some point, fluctuations become so strong 
that full turns of the spins are occurring regularly, 
making a description by a localized Gaussian un- 
suitable. This is reflected by the fact that the free 
energy has no absolute minimum but only a rela- 
tive one which vanishes at a critical temperature. 

We interpret this behavior as a phase transition 
by ergodicity: If we imagine each spin attached to 
a rubber band that accumulates the total rotation 
angle, fluctuations in the low-temperature phase 
do not reach over the crest of the cosine poten- 
tial, and spins will only rarely accumulate large 
rotations. The system is thus localized in a small 
region of phase space. In the high-temperature 
phase, spins frequently turn and will eventually 
reach any amount of rotation, spreading out the 
system periodically over the whole phase space. 

Wavelets do not describe the absolute direction 
of the spins and thus keep global 0(1) symmetry. 
They make it possible to obtain nonlocal correla- 
tions using only a small number of variational pa- 
rameters. However, the XY model is dominated 
by spin waves which are quite well described by 
plane waves, and the ability of wavelets to de- 
scribe complex irregular systems does not come 
into play in this problem. 

We find a critical inverse temperature of /3crit ~ 
0.69 for two dimensions on a 32 x 32-lattice and 
PcTit ~ 0.45 . . .0.49 in three dimensions. (The ac- 
curacy is limited by the shallowness of the relative 
minimum.) Low-temperature analysis of vortices 
predicts in two dimensions /9crit = S/tt 0.64. 

Fig. 1 shows the internal energy calculated by 
a Monte Carlo simulation and by the variational 
approach. 

4. Conclusions 

We have shown how wavelets can be used to 
describe fluctuations in a statistical field theory. 
Using a rather small number of variational pa- 
rameters, an approximate description of fluctu- 
ations in the XY model was obtained which is 
sufficiently accurate to predict a phase transition. 



3 1 




1 10 
beta 



Figure 1. Internal energy per site in the XY 
model vs. inverse temperature from variational 
wavelets (solid line) and from Monte Carlo calcu- 
lation (dotted line). 



The resulting phase transition temperatures are 
in good agreement with their expected values. 

Clearly, the approach is limited by the varia- 
tional principle used and the choice to expand in 
angle variables instead of spins. Further investi- 
gations are in progress to extend the treatment 
to SO{N) symmetric Landau-Ginzburg theories. 
As these theories are more easily accessible to 
position-space renormalization group transforma- 
tions, we are trying to represent such transforma- 
tions in wavelet space, effectively integrating out 
wavelet coefficients at the lowest scale. 
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